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THE ROLE OF STOCHASTIC ACCELERATION IN THE PROMPT EMISSION OF GAMMA-RAY BURSTS: 

APPLICATION TO HADRONIC INJECTION 

KOHTA MURASe'--, KATSUAKI ASANO", TOSHIO TERASAWA--\ AND PETER MESZAROs" 
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ABSTRACT 

We study effects of particle re-acceleration (or heating) in the post-shock region via magnetohydrody- 
namic/plasma turbulence, in the context of a mixed hadronic-leptonic model for the prompt emission of gamma- 
ray bursts (GRBs), using both analytical and numerical methods. We show that stochastically accelerated (or 
heated) leptons, which are injected via pp and 777 reactions and subsequent pair cascades, are plausibly able 
to reproduce the Band function spectra with a ^ I and /3 ^ 2-3 in the ^ MeV range. An additional hard 
component coming from the proton-induced cascade emission is simultaneously expected, which is compati- 
ble with observed extra power-law spectra far above the MeV range. We also discuss the specific implications 
of hadronic models for ongoing high-energy neutrino observations. 
Subject headings: gamma-ray burst: general — gamma rays: theory — radiation mechanisms: non-thermal 



1. INTRODUCTION 

The origin of prompt emission from gamma-ray bursts 
(GRBs) has been an open issue for more than 40 years. It 
is mainly observed in the E ^ E^'^ ^ 10 keV- 1 MeV range, 
most spectra being fitted by a smoothed broken power-law 
(the so-called Band function) with break energy £''" . The typ- 
ical low-energy photon index is a ^ 1 (where dF /dE oc £"" 
below ^ E^'^), while the typical high-energy photon index is 
/3 ^ 2-3 (where dF /dE oc £'~'^ above ^ £''"'). Observed light 
curves are highly variable. The variability can sometimes be 
^ ms, but pulses with various widths are observed. For typ- 
ical long GRBs, light curves consist of many variable pulses, 
so that the duration, T ~' 10'~ s, is usually longer So far, 
many models have been proposed to explain this most lumi- 
nous phenomenon in our un iverse (see recent reviews, e.g., 
iMesza-os 20061 IZhang 20071) . 

The classical scenario of the prompt emission is 
the optically-t hin sync hrotron internal shock model (e.g., 
iRees & Meszaros 19941) . where shocks are responsible for 
dissipation of the outflow kinetic energy and observed 
gamma rays are attributed to synchrotron radiation from 
relativistic electrons. However, there are a couple of 
issues in the classical scenario. First, the low-energy 
photon index may be incompati ble with values predicted 
by synchrotron emiss ion (e.g., iGhisellini & Celotti 19991: 
iMeszaros & Rees 20001) . This is the case especially when rel- 
ativistic electrons are in the fast-cooling regime, which leads 
to Q ^ L5, although the Klein-Nishina effect on synchrotron 
self-Compton (SSC) process may rep r oduce g^ 1 (e.g ., 
iDerishev et al. 2OOII: IWang et al. 20091 iDaigneetal. 201 ll) . 
One of the possible resolutions of this fast-cooling problem is 
to consider the magnetic field decay tim escale that is compa - 
rable to the electron cooling timescale (Pe'e r & Zhang 2006h . 
Injection of electrons to acceleration processes is one of 
the related issues. Shock acceleration of primary electrons 

' Dept. of Physics; Center for Cosmology and AstroParticle Physics, Ohio 
State University, 191 West Woodruff Avenue, Columbus, OH 43210, USA 

- Dept. of Physics. Tokyo Institute of Technology, 2-12-1 Ookayama, 
Meguro-ku, Tokyo 152-8550, Japan 

' Institute for Cosmic Ray Research, University of Tokyo, Kashiwa-no-ha 
5-1-5, Kashiwa-shi, Chiba 277-8582, Japan 

^ Dept. of Astronomy & Astrophysics; Dept. of Physics; Center for Parti- 
cle Astrophysics, Pennsylvania State University, University Park, PA 16802, 
USA 



may be inefficient for sufficiently small emission radii be- 
cause of the rapid cooling in the strong magnetic field r egion 
at the vicinity of the shock (iMedvedev & Spitkovskv 2009,) . 
If the outflow is magnetized enough, even for sub-luminal 
shocks, recent ab intio particle-in-cell simulations of col- 
lisionless shocks indicate that acceleration of nonthermal 
electrons at relativis tic shocks are inefficient whereas ions 
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can be accelerated (.Sironi & Spitkovskv 20 
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sue is the radiative efficiency problem (e.g.. lIoka et al. 20061: 
'Zhan g et al. 2007 ^. although this depends on estimates of the 
afterglow kinetic energy (Eichler & Waxman 2005). The high 
efficiency can be achieved by a large dispersio n in the Lorentz 
factor of the outflows (e.g.. .Beloborodov 200Q) . but it does not 
seem so easy t o reconcile that with t he observed spectral cor- 
relations (e.g.. lYonetoku et al. 20041) . A recent discussion of 
the classical scenario is found in Daigne et al. (201 1). 

On the other hand, many alternative models have been 
suggested. Diffusive synchrotron radiation can explain 
the low-energ y photon index if the fast-cooling prob- 
lem is fixed (iMedvedev 2000al) . Instead of dissipation 
of the kinetic energy via internal shocks, dissipation of 
the magnetic energy via, e.g., magnetic reconnection 
may play a crucial role as the dissipation mechanism, 
and observed gamma rays may be synchrotron or other 
electromagnetic radiation fro m rela tivistic elec trons (e.g. , 
|Spruit et al. 2001; Lvutikov_2QQ6t Zhang & Yan 201 It 
'McKinney & U zdenskv 20111) . If significant thermal en- 
ergy is stored in the GRB jet by energy deposition from 
the central engine and/or via some dissipation below the 
photosphere, the prompt emission may originate from quasi- 
thermal emission caused by Comptonization or Coulomb 



heating (e.g., [T hompson 1994; Ghiseflini & Celotti 1999 
Meszaros & Rees 2000: Pe'er etal. 200 6: loka et al. 2007 
Beloborod ov 20101) . One of the appealing possibilities is 
to consider the slow heating scenario. Around the photo- 
sphere, electrons that are slowly heated can lead to a ^ \ 
via Comptonization ([Thompson 19941: iPe'er et al. 2006h . 
where the typical electron energy is deter- 
mine d by the balance bet we en heating and cool - 
ing (" GhiseUini & Celotti 1999t IVurm & Poutanen 2009h . 
Alternatively, in the neutron-rich outflow. Coulomb col- 
lisions should work as slow heating below the baryonic 
photosphere, leading to hard photon spectra of a ^ -0.5 - 1 
and /3 ^ 2.5 via the pair injection by the np reaction and 
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subsequent cascades dBeloborodov 2010t IVurm et al. 201 lb . 
Such slow heating may be caused by stochastic accel- 
eration (so-called second-order Fermi acceleration) over 
the sub-hydrodynamical timescale and/or hydrodynamical 
timescale, which may occur in the magnetohydrody- 
namic (MHD) and plasma turbulences generated by 
shock or magnetic dissipation s uc h as magnetic reconnec- 
tions ("Bykov & Meszaros 1996t lAsano & Terasawa 20091 : 
IZhang & Yan 2011^ 

Emission at higher energies above ^ 10 MeV 
should also be important to reveal the radiation 
mechanism of the prompt emission, and various 
high- energy processes such as Comptonized thermal, 
SSC dRees & Meszaros 1994 lAsan o & Inoue 20 03), proton - 
induced cascade ( IVietri 1997: IDermer & Atovan 20061) . 
and proton synchrotron (iTotani 19981) emissions, were 
discussed. Such high-energy emission had been sparsely 
detected by the E GRET onboard the Compton Gamma-Ray 
Observatory (e.g., iHurleyetal. 19941) . Recently, however, 
significant observational progress has been achieved by 
Fermi. Dozens of GRBs have been observed by LAT onboard 
Fermi, and some bursts can be fitted _by the Ba nd function 
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up to GeV energ ies (lAbdo et al. 2009a 
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up to uev energ ies (lAbao et al. zuu^at lApgo et al. zuu^a 
lAbdo et al. 20101) . while some of them (GRB 090510, 
090902B, 0909 26A) clear ly have an extra hard component 

at > 10 MeV (lAbdo et a l. 2009b; Ackermann et al. 2010; 

— — ^ — 

ckermann et al. 201 Ih . GRBs with the extra GeV com- 
ponent belong to the brightest class of GRBs, which may 
indicate that this could be more com mon, even tho ugh it can 
be seen only in bright LAT GRBs (lGranot2010h . Another 
feature found by Fermi is that > 100 MeV gamma rays are 
delayed by ^ 0.1 - 1 s in the cosmological rest frame, behind 
the onset of the MeV emission. 

The origin of the high-energy emission is now under ac- 
tive debate. Late-time high-energy gamma-ray emission from 
GRBs such as 0809 16C, 090510, and 090902B has been at- 
tribute d to afterglow emission rather t han the prompt emis- 
sion jK umar & BarniolDuran 2OI0I; iGhisellini et al. 20101: 
iRazzaq ue 2011). But, the high-energy emission in the early 
phase usually has a strong variability and often correlates 
with the MeV em ission, which indicates an internal ori- 
gin jHe et al. 201 lb . Since the simple one-zone leptonic 
model often has difficulties, m u lti-zone models, no t only the 
SSC model (ICorsietal. 2010t iDaigneetal. 201 Ih but also 
the external inve rse-Compton (IC) model JToma et al. 2011*: 
iPe'er et al. 201 Ih and the synchrotron model (floka2010), 
have been invoked. The one-zone hadronic models are 
also viable, and the high-energy emission can be ex- 
plained b y proton-induced c a scade or proton sync hrotron 
emission (l Asano et al. 2009bt IRazzaque et al. 20101) . The 
hadronic models are of interest, since GRBs may be the 
sources of observed ultra-high-energy cosmic rays (UHE- 
CRs), whose origin h as been a big mystery for abo ut 50 
years (IWaxman 1995t IVietri 1995t iMurase et al. 20081) . and 
their neutrino signals may be seen by ongo ing neutrino obser- 
vations (a recent review is lWaxman 201 lb . 

In this paper, we study a slow heating scenario in the pres- 
ence of electromagnetic cascades initiated by injection at high 
energies. To explain the MeV emission self-consistenly, we 
consider effects of acceleration (via the second-order Fermi 
mechanism) or heating of secondaries, which can be antici- 
pated in the post-shock region of relativistic shocks (or mag- 
netic reconnection). We show that cascades can play an im- 
portant role in the injection of re-accelerated/heated particles. 



and slow heating scenarios via stochastic acceleration or wave 
heating typically lead to hard spectra compatible with obser- 
vations. As an injection process at high energies, we consider 
an application to hadronic injection via pp and pj reactions, 
whereby the hadronic models have the appealing feature of 
explaining an extra hard component at GeV energies. 

The paper is organized as follows. Section 2 describes an 
overview of the model, and we provide the numerical method 
and its results in Section 3. We discuss implications for neu- 
trinos in Sections 4, and our results are summarized in Sec- 
tion 6. Throughout this work, cosmological parameters are 
set to //() = 71 km s"' Mpc~', Qm = 0.3, and J7a = 0.7, and we 
adopt the conventional notation Q = Q^ x 10'. 

2. THEORETICAL MODEL 

In this work, we investigate the role of particle re- 
acceleration, motivated by the goal to resolve the low-energy 
photon index and fast-cooling problems. In the post-shock 
region of relativistic shocks, turbulent magnetic fields would 
be generated, where stochastic acceleration (or second-Fermi 
acceleration) may work efficiently since the Alfven velocity 
is expected to be close to ^ (0.1 - l)c. The second-order 
Fermi acceleration over the sub-hydrodynamical scale plays 
a role of slow heating to avoid the fast cooling problem 
of relativistic electrons, and it may give harder electron 
spectra than the typical first-order Fermi acceleration. This 
possibility was proposed by Bykov & Meszaros (1996) 
and recently studied by Asano & Terasawa (2010). Also, 
particles may be slowly heated by MHD waves, as has 
been discussed in the solar co rona (e.g.. [ Alf ven 1947 

'Fossum & Carlsson 2005"; jSuzuki & Inutsuka 2005 

Tomczyk et al. 2007; Jess et al. 20091) 

Such slow heating operates not only for primary electrons 
but also electron-positron pairs generated via cascades initi- 
ated by some high-energy injection. Motivated by the elec- 
tron injection/acceleration problem and the existence of an 
extra hard component at > 10 MeV, this work focuses on 
a hadronic model, where pairs are produced via the electro- 
magnetic cascade initiated by hadronic injection from pp(np) 
and /?7 reactions, i.e., the proton-induced cascade. Although 
there is a similarity to the collisional scenario proposed by 
Beloborodov (2010), our scenarios are different in that lep- 
tons are heated by turbulence rather than Coulomb heating 
and the dissipation radius can be much larger. 

The hadronic emission associated with prompt emission 
was originally motivated by t he hypothesis t hat observed 
UHECRs come from GRBs JWaxman 19951; IVietri 19951; 
iMurase et al. 20081) . This requires that the UHECR en- 
ergy (at ~ 10 eV) is comparable to the gamma-ray en- 
ergy, £hecr = E\^{dNpldE\,) - £^. Although the GRB 
rate evolution, total gamma-ray energy and proton spec- 
tral index are still uncertain, the total nonthermal baryon 
loading i s expected to be large, £cr = ^-cs.l^i ^ 3 - 
300 (e.g.. IMurase et al. 2008t IWaxman 20101 and references 
therein). If sufficiently large baryon loading is realized 
and the photomeson production efficiency is high enough, 
the hadronic emission can dominate over the leptonic emis- 
sion at high energies (lAsano et al. 2009al) . and the visibil- 
ity of the hadronic components at high energies is enhanced 
when the high-energy photon spectrum is much steeper than 
/? ^ 2. The extra hard component in the GeV range, which 
was s een by Fermi, can b e explained by the hadronic emis- 
sion ( lAsano et al. 2009bl) . But, in the previous work, the 



The Role of Stochastic Acceleration in the GRB Prompt Emission 



Band function was assumed ad hoc, otherwise the low-energy 
photon index becomes too soft (which leads to overestimating 
gamma rays and neutrinos). 

In the next subsection, we consider shock dissipation in 
the baryonic outflow and shock acceleration of particles. 
Then, we discuss stochastic acceleration (or heating) by 
plasma/MHD turbulence in the downstream region to take 
into account its effects with a simplified model. 

2.1. Shock Acceleration at Internal Shocks 

In the classical scenario, electrons are accelerated at in- 
ternal shocks, and the observed gamma-ray emission is at- 
tributed to electromagnetic radiation from non-thermal elec- 
trons. If the outflow contains baryons, it is natural to expect 
that protons are also accelerated on a timescale 

facc = ?7— , (1) 

C 

where r? ^ 1 - 10 in th e efficient case (e.; 
iBednarz & Ostrowski 19961: iRachen & Meszaros 1998 
The acceleration timescale is typically fast, and pro- 
tons that are accelerated up to sufficiently high en- 
ergies should interact with non-thermal photons via 
photo meson production, le ad ing to ^ PeV neutri - 
nos (IWaxman & Bahcall 1997 ': 'Rache n & Meszaros 1998h . 
GeV- TeV gamma rays (IVie tri 1997), and UHE gamma 
rays (Murase 2009). Around the photosphere, GeV-TeV neu- 
trinos and GeV-TeV gamma rays are also produced t hrough 
the pp(np) reaction (iMurase 2008HWang & Dai 20091) . 

In our model, we consider a similar picture in which 
charged particles are accelerated at shocks. Motivated by the 
electron injection/acceleration issue, however, we focus on 
the case of hadronic injection (neglecting primary electrons), 
where proton-induced gamma-ray emission is especially rel- 
evant. When the shock dissipation radius is small enough, the 
efficient ppinp) reaction occurs in the baryon-rich outflow, 
and high-energy gamma rays and pairs are produced as well 
as neutrinos. Also, sufficiently high-energy protons can inter- 
act with soft photons which may co me from residual therma l 
emission around the photosphere (lAsano & Takahara 2003h 
(and/or non-thermal emission from primary electrons if their 
acceleration occurs). Those pp and 777 reactions lead to 
production of high-energy gamma rays with E ^ Q.lEp ~ 
TeV £^,12 and pairs with E^ ^ Q.QSEp ~ 0.5 TeV Ep\2. 
Lower-energy photons should also be produced via the syn- 
chrotron and IC processes, which eventually affect the distri- 
bution of seed photons so that the processes are nonlinear. 

Sufficiently high-energy gamma rays should lead to pair 
creation with seed photons at EE « T^mlc*, so that the fol- 
lowing cascade is unavoidable. As long as the pair-creation 
opacity of high-energy gamma rays is larger than unity, 
gamma rays lead to generation of pairs that lose their ener- 
gies via synchrotron and IC emissions. As a result, the cas- 
cade gamma-ray spectrum will be eventually formed, where 
the number of pairs can be larger than that of pairs injected 
directly from the pp and pj reactions. 

We are interested in cases where the pp or 777 reaction 
efficiency is high enough, otherwise the hadronic compo- 
nent is insufficient to explain the high-energy gamma-ray 
emission observed by Fermi. Hence, for the purpose of 
demonstrations, this work focuses on relatively small dissipa- 
tion radii of r ^ 10'"~' cm such that tj- « arnthl ^0.1-1 
(where nth is the thermal proton density and / is the co- 
moving length), though large uncertainty allows us to con- 



sider different cases wher e the pp and P7 re action efficien- 
cies are much lower (e.g.. lMurase et al. 2008h . The pp reac- 
tion effi ciency is rough ly estimated to be fpp « Kppappn^hl — 
O.OStt- (iMurase 20081) . where Kpp is the pp inelasticity of 
protons. We expect that the observed gamma-ray spectrum 
will be eventually formed (see the next section). Then, as- 
suming the (broken) power-law photon spectrum, the ef- 
fec tive photomeson product i on efficiency is estimated to 
be (IWaxman & Bahcall l"997HMurase & Nagataki 2006bh 

^7.51.5 f (Ep/E^;f-' (Ep < E);) 

n3ri,(£>'7500keV) \ (Ep/E^^f-' (E^; < Ep). 

(2) 
where E^' ~ 80 PeV (£''7500 keV) ^rjj is the energy of pro- 
tons interact with photons with E^'^ and the multi-pion produc- 
tion effect is included. 

Gamma rays also interact with the same target photon 
field, so that the above equation implies that the pair-creation 
opacity is (e.g.. 'Lith wick & Sari 200!t IMurase & loka 200"8l: 
l&ipta & Zhang 2008) 
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where E'^' ~ 130 GeV (£fa7500keV) Tj-j and f(J3) - 0.1 
for /3 ~ 2 (IBarin g 2006). (Note that the opacity in the high- 
est energies can be smaller than unity due to existence of the 
synchrotron self-absorption cutoff.) Then, in a simple one- 
zone model, the pair-creation break (or cutoff) is estimated 



to be £■'= 



2/3 



2-li 
/3-1 



Most of the high-energy gamma rays are radiated below this 
energy, and the above equation suggests that the effective in- 
jection rate of electron-positron pairs is 



eg, 



min[l,</ines >]UpGe 



t& 



(4) 



yn 



where Up is the total cosmic-ray proton energy density, fdyn ~ 
3//c is the dynamical time, < /mes > is the effective meson 
production efficiency averaged over proton energies, and G^ 
such that J deO^ = 1 is determined by the details of the pair 
cascade. Because of the cascade and the multiplicity of the 
hadronic reactions, the numb er of pairs can be even larger 
than that of the reactions (cf. Beloborodov 20101) . If there is 
no re-acceleration/heating and escape, we get 



dn^ 



d_ 
de 



de 
— I 
dt 



+ Qh 



where 



de A 2 

-r = -^<^tc(Ub+Ukn)(i 
dt 3 



1) 



(5) 



(6) 



and Ub is the magnetic field energy density and Ukn is 
the photon energy density with the correction by the Klein- 
Nishina effect. Therefore, at energies where injection occurs, 
we expect that the cascade typically leads to a flat energy 
spectrum of electrons, so that we obtain EF{E) ex £'0~''-5 that 
is also seen by our numerical calculations. Such hard photon 
spectra can explain the extra hard compone nt at > 10 MeV 
observed by Fermi (e.g.. lAsano et al. 2009bl) . In the hadronic 
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model that attributes the extra high-energy component to the 
proton-induced cascade, one typically expects 



'min[l,</mes >]fcR, 



(7) 



where £^ex is the energy in gamma-rays of the extra compo- 
nent, which is typically smaller than the total gamma-ray en- 
ergy released as prompt emission, 8^. As the emission radius 
is larger (leading to lower < /mes >), larger £cr is required. 

2.2. Slow Heating by Turbulence 

Next, we consider stochastic acceleration of particles that 
are produced by injection at high energies and subsequent 
cascades. The turbulent magnetic field in the downstream 
of the shock is expected to be significantly amplified. Cas- 
cade pairs that are rapidly cooling (i.e., they are in the 
fast cooling regime) are distributed in the downstream re- 
gion, and can be re-accelerated by plasma/MHD turbu- 
lence via the second-order Fermi mechanism. Assuming 
isotropy of the particle momenta and of the fluctuations, one 
may apply the Fokker-Pla nck equation (e.g.. iLiu et al. 2006t 
iBrunetti & Lazarian 20071) . 



dn 



dn. 



— \ D 

dt de \ "" de 



d_ 
de 
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■7^ + 2/, (8) 
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where A = (D£e/£)(1 -7-2)/(l +7"') and D^e is the dif- 
fusion coefficient. For isotropic magnetic field fluctua- 
tions with spectral energy density Wu oc k~^, one may write 
Dee oc e^ with q = q for resonant acceleration, where scat- 
tering and acceleration timescales of charged partic les are 
proportional to e^I Dee oc e^"'' (e.g.. Miller & Ramatv 19891: 
iDermeretal. 1996|). For non-resonant accele ration, one may 
expect Dee oc e^ (IBrunetti & Lazarian 2007h . The fluid and 
particles behind the shock are advected to the downstream. 
The length scale of the turbulent region will be Ztui ~ (c/3)ftiir, 
where ftui is the typical lifetime of the turbulence. Then, since 
the Alfven velocity can be close to ~ (0.1 - l)c given that 
the magnetic field is significantly amplified (/34 = B/^/A-kIi ^ 

^(3Ub/2U±) — 0.4eB-i, where h is the enthalpy density 
and eg is the energy fraction carried by magnetic fields), 
for fluctuations with a simple power law, the stochastic 
acceleration timescale is expressed as (jDermer et al. 19961 : 
iBykov & Meszaros 1996HPetrosian & Liu 2004 



4De 



— Jysta " 



2-q 



(9) 



where r/, is the Larmor radius, rysta ~ 1 is a pre-factor that 
depends on the wave amplitude and the Alfven velocity, 
and we have assumed that In/k^in ^ Im- This stochastic 
acceleration timescale can be longer than the shock accel- 
eration timescale, face = Vl/c- But it can still be shorter 
than the hydrodynamical timescale, /mr/c, in the relevant 
energy range where ri ^ /tur (except for exactly q = 2 
over all k). In this sense, the stochastic acceleration plays 
the role of a slow heating over the sub-hydrodynamical 
timescale ( Bvkov & Meszaros 19961) . 

Various numerical simulations, includ ing plasma sim- 
ulations for colhsi onless shocks (e.g ., ISpitkovsky 20081: 
IChang et al. 20081: iKeshet et al. 20091) and MHD sim- 
ulations (e.g., |3iang et al. 2003 iMizuno et al. 201 It 
llnoue et al. 201 1[) . have recently been attempted. For the 
MHD scale turbulence, Inoue et al. (201 1) showed that highly 
relativistic turbulence decays fast, whereas the transonic and 



subsonic turbulences are maintained much longer, and 
their timescales are determined by the initial scale of the 
inhomogeneity in the fluid density, which may be comparable 
to the shell width. Such long-lasting large scale turbulences 
may play an important role to accelerate/heat high-energy 
electrons (and positrons) that emit MeV gamma-rays. 

On the other hand, small scale turbulences induced via 
plasma instabilities, such as seen in PIC simulations (e.g., 
ISpitkovskv 20081: IChang et al. 20081) seem to settle quickly 
and exist only in the vicinity of the shock at earlier times. 
Keshet et al. (2009) suggested high-energy particles play an 
important role for the long term evolution and the magnetic 
field is progressively generated on larger scales in a growing 
region around the shock, where colhsionless shock configu- 
rations do reach the steady state. Effects of high-energy par- 
ticles may be important for the generation and sustainment 
of magnetic field turbulences. The growth of the turbulence 
scale, as seen in coalescence of current beams for the case of 
the Weibel instability (Kato 2005), should be tested through 
long-term simulations including the non-thermal particles. 

Below, we mainly expect that the resonant acceleration 
via scatterings with Alfven waves can work, whereas one 
might expect other possibilities such as the non-resonant 
acceleration or heating by plasma/MHD waves. However, 
the non-resonant acceleration due to compressive modes 
will not work for (3^ > /3cs (sound speed divided by c) 
since the fast and slow modes will be damped due to 
the Landau damping, although one may have possibili- 
ties of some non-resonant scattering by Alfven waves or 
magnetic bottles formed by large-scale slow-mode per- 
turbations for (3a < /3cs ~ l/\/3 (due to small values of 
e^) dBeresnvak et al. 201 lb ^. On the other hand, dissipation 
of MHD waves produced at the shock can play crucial roles 
on heating particles. Heating of par ticles has been debated in 
the context of the sol ar c orona (e.g..lSuzuki & Inutsuka 20051 : 
iTomczvk et al. 20071 iJess et al. 20091) and a dvection- 
dominated accretion flow (e.g.. lMedvedev2000bh . but the 
details are uncertain in the case of the relativistic plasma 
expected for the GRB shock. Though the detailed study is 
beyond the scope of this work, we expect the analogous slow 
heating effect on forming a spectrum as long as particles in 
the relevant energy range can be heated by the waves. 

At present, we still have no complete model of the tur- 
bulence for GRB internal shocks. Hence, as in previ- 
ous work (Asano & Terasawa2009), we consider the sim- 
plest setup where the stochastic acceleration time is energy- 
independent (Dee oc £^) in the range in which we are in- 
terested. Such a situation may be realized for the reso- 
nant acceleration, if g' ~ 1.5 — 2 at relevant electron ener- 
gies such that ri ^ /tur, and necessary steep turbulent spectra 
can be expected for the MHD t urbulence (q ^ 1.5-2) (e.g., 
ICho et al. 2003tlLiu et al. 20^6h or c oUisionless shock s in the 
presence of cosmic rays (q ^ 2) (e.g.. lKatz et al. 2007h . 

Note that the turbulent spectrum may not be a simple 
power law so q in the diffusion coefficient is rather energy- 
dependent in the wide energy range, and the MHD simula- 
tions have suggested the spectrum is flat at large scales but 
steep enough at small scales above the transition scale (e.g., 

^ The Landau damping rate would be small if the distribution function 
of thermal ions have negligible slopes, df/dv, around the resonant velocity. 
But, in this case, the majority of ions could be trapped in finite amplitude 
slow-mode waves, so possible nonlinear effects should be incorporated to 
obtain the damping rate. 
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IChang et al. 20081: llnoueetal. 2011b . But it would be ap- 
proximately valid to assume a singe value of q in the relatively 
narrow energy range relevant for the MeV emission. Al- 
though future numerical studies on plasma/MHD turbulence 
are required for more realistic calculations, our phenomeno- 
logical approach is enough to demonstrate the role of stochas- 
tic acceleration and we show that the process efficiency is rea- 
sonably sufficient to explain observations given the setup. 

Next, we briefly discuss the expected spectra of particles 
re-accelerated via stochastic acceleration in the plasma/MHD 
turbulence, although the details depend on diffusion, en- 
ergy losses, and particle escape. Let us consider low 
energies where energy losses are negligible. Then, for 
q = 2, the steady-state spectrum is obtained as n^ oc 
g0.5-(9/4+4/,„/re.c)'/- jLiu et al. 20061) . At sufficiently low ener- 
gies (ri ^ /tui) that are relevant in the case of GRB prompt 
emission, fgsc is so long that we may drop it, where one has 
«£ ex e'"'. This means that the stochastic acceleration leads to 
hard electron spectra, and the synchrotron spectrum becomes 

EF(E) ex E^-^/^ (E < E^'), (10) 

for q > 4/3 (and we get EF(E) ex £"^/^ for q < 4/3). Hence, 
the low-energy photon index of a ~ q/2 ^ 1 is possible if ^ ^ 
1.5-2 (which seems typical for the MHD turbulence). In the 
classical scenario of prompt emission, electrons are typically 
in the fast cooling regime, so that high-energy electrons cool 
down within the advection time, ^ 3l/c. But, this problem 
can be avoided when acceleration (or wave heating) plays a 
role of slow heating. 

When cooling is absent, the maximum energy of re- 
accelerated/heated electrons should eventually be limited, 
since they can carry only a fraction of the turbulent energy. Or, 
radiative cooling is so strong that the typical electron energy is 
determined by the balance between acceleration/heating and 
cooling. If we knew the turbulent spectrum from first prin- 
ciples, we would be able to calculate 7typ for the stochas- 
tic acceleration. For example, if the turbulent spectrum is 
flat at large scales and then becomes sufficiently steep at 
^ lOOfemin, the Stochastic acceleration time can be fsta ~ 1 .3 x 
10"'' s (/tui/0.20'"i3r2 7. Comparison with the synchrotron 
cooHng time of fsyn =i 1.8 x 10"^ s eB]_ir^i^rljL;^^^^^'y2j (ith 
is the luminosity associated with the thermal energy) gives 
7typ - 690 e-;_,(Ztur/0.2/)-Vi3ri7L7,ig3 5 leading to ^ MeV 
synchrotron gamma rays. Unfortunately, however, the realis- 
tic stochastic acceleration rate is highly uncertain at present so 
that we instead make an ad hoc assumption that the stochastic 
acceleration rate is determined by requiring that the observed 
break energy is the synchrotron energy of electrons with 
the typical electron Lorentz factor (lAsano & Terasaw a 2009,), 
where we find*" 



7typ' 



cbi- \ 1/2 



TheB 



1/4 f_ 



■^^""i^'^-'lsOOkeV 



1/2 



-1/4 
-"th,53.5- 



(11) 



Phenomenologically, one can obtain another expression, 7typ 



for <? < 2 or 7,yp ~ ln(:2^)| ^7p,h for g ~ 2. 



3-g €, mp l-q 

Here e^ is the energy fraction carried by electrons accelerated via turbulence 
and 7icacc is the electron Lorentz factor at which the stochastic acceleration 
does not operate (which can be expected when particles do not interact with 
plasma/MHD waves). Then, fe is the effective number fraction of electrons 
that are injected to the stochastic acceleration, which is determined by the 
details of the cascade and acceleration processes. 



As described in the next section, we take into account slow 
heating (stochastic acceleration or wave heating) of leptons 
by using the Monte Carlo method, assuming that the accel- 
eration/heating time is energy-independent and equation (11) 
applies. Although it is a toy model, we can demonstrate that 
the re-accelerated/heated particles have hard spectra leading 
to a ^ 1 , as expected in this subsection. 

3. NUMERICAL RESULTS OF GAMMA-RAY SPECTRA 

In the previous section, we gave the basic picture of our 
model with analytical expressions. But many microphysical 
processes have to be taken into account in order to calculate 
gamma-ray spectra. O ne of the appro aches is to solve the 
kinetic equations (e.g., ISvensson 1987h . but it is time con- 
suming for our problem. Alternatively, we employ the Monte 
Carlo code used in, e.g., Asano et al. (2009a), Asano & Tera- 
sawa (2009), where the energy distributions of all particles are 
simulated iteratively until they converge to a self-consistent 
steady state, which is assumed to be realized within the pulse 
timescale. 

All the important microphysical processes for the electro- 
magnetic cascade, 77 pair creation by electron and positron 
pairs, synchrotron, and Compton (both in the Thomson 
and Klein-Nishina regimes) emissions from all relativistic 
particles, an d synchrotron self-abso rption, are properly in- 
cluded (e.g., iBlumenthal & Gould 1 970). The photomeson 
production is approximately included as in Asano & Na- 
gataki (2006), using experimental data but neglecting high- 
multiplicity/inelasticity, multi-pion production at high ener- 
gies, and it is assumed that a muon neutrino from pion decay 
carries a quarter of the pion energy, and neutrinos and a lepton 
from muon decay carry a third of of the muon energy. Note 
that such a simplified treatment is timesaving and sufficient 
for this work, since effects of more detailed microphysics are 
typically moderate dMurase & Nagataki 2006al:lMurase 20071: 
iBaerwald et al. 201 Ih . The Bethe-Heitler process by protons 
is also included using the cross section and inelasticity given 
by Chodorwski et al. (1992), but this effect is only modest 
in the re levant energy rang e for GRB prompt emission spec- 
tra (e.g.. iMannheim et al. 2000 ). In addition, we implement 
in detail the pp reaction for which high multiplicity and in- 
elasticity in pion production is important. The high-energy 
process is calculated as in Murase (2008), whereas the low- 
energy p rocess below 100 GeV is included based on runs of 
Geant 4 jAgostinelli et al. 2003h . The Thomson scattering ef- 
fect due to thermal electrons, which would weaken photons of 
energy above ~ irigC^ in the comoving frame, is included by 
assuming that photons interact with electrons of temperature 
100 keV within a timescale l/c before they escape from the 
emission region. The spectral softening by this effect has been 
calculated with a Monte-Carlo method with the Klein-Nishina 
cross section calculated in advance. 

There are several input parameters required for calcula- 
tions. The emission radius r and the bulk Lorentz fac- 
tor r are relevant parameters. One typically expects r ~ 
lO' - lO'^^ cm in the classical internal shock scenario, or 
r ^ lO" - lO' cm in the baryonic photospheric scenario. We 
consider relatively small radii (but larger than the coasting ra- 

12 1 ^ 

dius), r ~ 10 - 1 cm, wher e the efficient hadronic reac - 
tions are efficient (' Asano 20051 : iMurase & Nagataki 2006ah . 
We also take F ^ 100- 1000, and the comoving length is given 
by I = r/T. Another important parameter is the luminosity of 
protons accelerated at the shock, L^, = AnP'T^cUp (and the 
nonthermal baryon loading parameter is expressed as ^cr = 
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Fig. 1. — Gamma-ray spectra of proton-induced cascade emission without 
stocliastic re-acceleration. The source redshift is assumed to be z = 1 . The 
strong magnetic field is assumed over the length scale of /tm =0.17/ and rj=\ 
is assumed. The corresponding Thomson depth is tt = 1 . 



Lp/L^). In this work, we consider Lp 



10 



53 



10^^ erg s-\ 



assuming the proton distribution with dnp/ds ex e. e ^^p'^p ' 
for Ep > 7th'«pC^ = 4m pC^. The maximum energy e™^" is 
given by the comparison between the acceleration time and 
cooling timescales as well as the Hillas condition. We also 
assume the luminosity associated with the thermal energy, 
Lth = 47rr^r^c[/th(^ Lp), where Uth = JihnihnipC^ is the thermal 
energy density. The comoving magnetic field is introduced as 
%(= efitAh) ~ 0.01 - lUp. In accord with our motivation and 
purposes, we neglect the acceleration of primary electrons. 

In Figure 1, we show an example of the proton-induced cas- 
cade emission. The pj and pp reactions lead to very high- 
energy gamma rays which eventually cascade down via pair 
creation, synchrotron and Compton emissions. (Note that we 
include both IC scattering and Compton downscattering.) The 
proton synchrotron component also contributes to the intrinsic 
very high-energy component but it is significantly attenuated 
and cascaded in the source. Then, as expected before, a flat 
and hard spectrum, which is different from the Band function, 
is typically obtained below the pair-creation cutoff. Note that 
such a flat and hard spectrum is expected for larger dissipation 
radii as studied in detail by Asano et al. (2009a). 

Now, we take into account effects of stochastic accelera- 
tion (or heating). We give the energy gain/loss per collision, 
Ae/e, following the Gaussian probability function, according 
to Asano & Terasawa (2009) that focused on the case of low- 
energy injection where the cascade does not play an important 
role. Then, we set the stochastic acceleration (or heating) rate 
such that the synchrotron peak of 7typ electrons correspond 
to £''"'. We also introduce sub-parameters, 7reacc and Ztur, but 
our main results are not sensitive to those additional parame- 
ters. The former is related to how long stochastic acceleration 
(or heating) operates, and 7reacc ^ 7typ is needed for the re- 
acceleration/heating to work well. Since the calculations for 
7reacc ^ 1 are too time-consuming, we take 7i-eacc ^10- 100, 
which is enough for our purpose. The latter parameter is 
physically motivated by the finite lifetime of the turbulent 
magnetic field, and we assume l^^i/l ^0.1-1 in this work. 
Thus, our "one-zone" code divides the shocked region into 



two parts: a strongly magnetized region of scale /tur and the 
second, remaining part of scale / - Ztur- The photon density is 
assumed to be homogeneous over the entire scale I, and sec- 
ondary particles from protons would be injected also homoge- 
neously in this scale. But, for simplicity, we treat only a frac- 
tion /tui// of the secondary particles injected in the strongly 
magnetized region. Although IC emission by the particles 
injecte d in the non-disturbed re gion may appear in the GeV 
range (lAsano & Terasawa 2009 ). we neglect such emission. 
The absence of heating in this region may make such emis- 
sion less important. The parameter /tur can affect high-energy 
IC spectra quantitatively, but the qualitative features of our re- 
sults are not altered (see below). The calculation time is set to 
?tur = 3/tur/c for the secondary particles, while that for acceler- 
ated protons is //c. We do not consider stochastic acceleration 
of protons since more efficient shock acceleration is assumed. 
We also remark that the energy ranges for pion-producing pro- 
tons and electrons with ^ 7typ are quite different. Also, we 
neglect re-acceleration of pions and muons, though it might 
be potentially relevant (see appendix A). 

In Figure 2, we show the case where stochastic re- 
acceleration (or heating) is included, for the same parame- 
ter set shown in Figure 1. When the re-acceleration/heating 
is turned on, electrons and positrons can avoid their fast 
cooling, and the energy balance between acceleration/heating 
and cooling gives the typical energy of 7typ ^ 100. Syn- 
chrotron emission forms another component in addition to 
the flat proton-induced cascade component, and its peak can 
be attributed to the observed peak energy. The non-linear 
effect on the photomeson production enhances the proton- 
induced cascade component compared with the case in Fig- 
ure 1. As argued in Section 2.2, one sees that a hard elec- 
tron distribution indeed leads to the low-energy photon in- 
dex of a ^ 1. The high-energy synchrotron spectrum of 
re-accelerated/heated particles is expected to become steeper 
and steeper above the synchrotron peak. But, as a result of 
the superposition of the slow heating component and the un- 
derlying cascade component, the high-energy photon index 
is effectively regarded as /3 ^ 2.5 at ^ 1-10 MeV ener- 
gies. As demonstrated in Figure 2, in the slow heating sce- 
nario with the electromagnetic cascade initiated by some in- 
jection at high energies, the resulting spectrum consists of 
an apparent broken power-law component due to the syn- 
chrotron emission from re-accelerated/heated particles and an 
extra hard component from the cascade emission (including 
IC emission by re-accelerated/ heated particles), which i s dif- 
ferent from the previous work (lAsano & Terasawa 20091) that 
instead obtained a simple power law with a high-energy cut- 
off. Such a spectrum may also be responsible for observe d 
Fermi GRBs, e.g., GRB 090926B (lAckermann et al. 201 Ih . 
Note that the gamma-ray luminosity in the MeV range comes 
from re-accelerated/heated leptons whose energies are sup- 
plied by the plasma/MHD turbulence, so that it can be larger 
than a fraction of the proton luminosity that is transferred to 
electromagnetic components via decay of mesons and muons. 
The spectral break at ^ 10 keV is determined by the parameter 
7reacc, bclow which the stochastic acceleration (or heating) be- 
comes ineffective. This parameter is kept as a free one in this 
simulation so that the quantitative value of this low-energy 
break does not carry strong implications. This break's ap- 
pearance is interestingly similar to the observed one reported 
in some of Fermi GRBs. 

In Figure 3, we show the fate of shock-accelerated high- 
energy protons that are injected with a simple power law. The 
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Fig. 2. — Gamma-ray spectra of proton-induced cascade emission with 
stochastic acceleration (or heating). The source redshift is assumed to be 
z = 1. Used sub-parameters (that are not critical for the result) are /tm// = 
0.17, 7rcacc = 10, and r? = 1. The obtained gamma-ray luminosity is L-, ~ 
1.2 X 10^^ erg s"' and the corresponding Thomson depth is tt = 1. 



photomeson production efficiency is so high in this param- 
eter set that almost all the protons are eventually depleted. 
Note that UHECR production wit hout their depletion is pos- 
sible at larger d issipation radii (IWaxman & Bahcall 19971: 
iMurase et al. 200 8) but larger baryon loading would be re- 
quired to achieve the same hadronic gamma-ray flux. One 
sees that high-energy mesons are mainly produced by the pho- 
tomeson production, while lower-energy mesons come from 
the pp reaction, and the result is in agreement with the an- 
alytical expectation described in Section 2.1. Hence, when 
protons are accelerated up to sufficiently high energies, one 
expects that the photomeson production is typically dominant 
in cascaded gamma rays, though the signature of pp reactions 
can be seen in TeV neutrinos for sufficiently small dissipation 
radii. However, the pp reaction can be more imp ortant \y hen 
there is no proton acceleration (e.g., Beloborodo v 20101) . or 
the nonthermal proton spectral index is steeper than 2 or pro- 
ton acceleration is inefficient (Murase 2008.) . Even in these 
cases, our model can work to reproduce the observed MeV 
emission. Such a case is demonstrated in Figure 4, where we 
obtain a broken power-law spectrum that can be compatible 
with the Band function. The value of rj used here (face.;; = 
rj{Ep/TeB)) seems much larger than conventional ones that 
are required for UHECR acceleration, although such values 
have b een use d to inte rpret observed blazar spectra (e.g., 
Ilnoue & Takaha ra 1996h . 

In Figure 5, we show the result for larger emission radii, 
taking into account effects of stochastic acceleration (or heat- 
ing). In this figure, a smaller magnetic field is assumed, so 
that IC emission from re-accelerated/heated leptons, which is 
expected around 5 GeV 7^^ 2(^'"/500 keV), become more 
prominent. The effect of the IC emission is generally non- 
linear, and the dependence on Ztur is demonstrated in Figure 6. 
One can see that, for larger values of /tur, more particles are 
re-accelerated/heated and lose their energies via the IC emis- 
sion. As a result, more gamma rays are attenuated and cas- 
caded, which enhance the extra hard component and also pro- 
vide more seed photons for the photomeson prediction. The 
different efficiency for IC emission results in the variety of 
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Fig. 3. — Spectra of injected and survived protons in the outflow comoving 
frame for the calculation in Figure 2. The resulting pion spectra are also 
shown. 
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Fig. 4. — Gamma-ray spectra for the case where proton acceleration is 
inefficient. The source redshift is assumed to be z = 1 . Used sub-parameters 
are hur/I = 0.99, 7rcacc = 10, and t] = lO''. The obtained gamma-ray luminosity 
is L-y ~ 1.6 X 10^^ erg s"' and the corresponding Thomson depth is T7- = 1. 



the high-energy spectral index /3, which is from 2.8 to ^2 
in Figure 6. One also sees that Ztur cannot be too small to 
have sufficiently large gamma-ray luminosities in the MeV 
range, but the result is not so sensitive as long as Ztur is large 
enough. Note that, the other sub-parameter, 7reacc^ affects the 
flux ratio of the broken power-law component to the power- 
law-like cascade component rather than /tui, since more lep- 
tons can be injected to the re-acceleration/heating process for 
smaller 7reacc- We do not show such other cases because of 
the onerous calculation time involved, but our qualitative re- 
sults are not much altered by these sub-parameters as long 
as 7reacc ^ 7typ and Ztur is large enough. This is because the 
timescale for remaining above 7reacc would be close to the dif- 
fusion timescale in energy space, ^ e^ /{AD^^) = fsta- There- 
fore, it would not be extremely large (E^gj^/E^ > 0.01), un- 
less the contribution from primary electrons become relevant. 
Although we have considered relatively small emission radii, 
such a spectrum composed of an apparent broken-power law 
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Fig. 5. — Gamma-ray spectra of proton-induced cascade emission with 
stochastic acceleration (or heating). The source redshift is assumed to be 
z = 1. Used sub-parameters are /tur// - 0.36, 7reacc = 100, and r? = 1. The 
obtained gamma-ray luminosity is L-y ~ 8.6 X 10^' erg s"' and the corre- 
sponding Thomson depth is t?- = 0. 1 . 
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Fig. 7. — The total neutrino fluences (for all flavors) from a ORB event 
at z = 1 . The dot-dashed, dashed and solid lines represent the fluences con'e- 
sponding Figures 4, 3, and 5, respectively, in order of the relative importance 
of the pp reaction. 
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Fig. 6. — Same as Figure 5, but dependence on the sub-parameter, /tm, is 
shown. 



and an extra hard component is expected at larger radii too 
as long as high-energy injection and efficient cascades are 
achieved. 

4. SPECIFIC IMPLICATIONS OF HADRONIC MODELS FOR 
HIGH-ENERGY NEUTRINO OBSERVATIONS 

The purpose of our paper is to investigate the role of 
stochastic acceleration (or wave heating), and we specifi- 
cally consider hadronic processes as high-energy injection 
processes. At the same time, it would be useful to discuss 
implications of hadronic models that can explain observed 
gamma rays. One of the important tests of the hadronic 
models is the prediction of high-energy neutrinos. The 
high-energy neutrino emission is a generic consequence of 
hadronic models which explain observed gamma-ray spectra 
through proton-induced cascade emission rather th an via pro- 
ton synchrotron emission (e.g., iBecker et al. 20 IQ) . whether 
the re-acceleration/heating occurs or not. In the hadronic 



model explaining the extra hard component with the proton- 
induced cascade, one typically expects, Ei, ^ Ejex- This rela- 
tion between E^, and ^..^ex is not much affected by the multi- 
pion production effect, even though estimates of /,nes and the 
required Ecr can be affected. However, we do not know typi- 
cal values of Eyex at present, and recent analyses indicate that 
the GeV c omponent is less than 1 % of the MeV one on av- 
erage (e.g.. lBeniamini et al. 201 Ih . 

In Figure 7, we show neutrino spectra for three param- 
eter sets shown in the previous section, when meson re- 
acceleration is irrelevant. For normalization, the total gamma- 
ray energy is set to El^° =10 '^ erg and z = 1 is assumed. For 
the thick curve where rj = 0.1, PeV-EeV neutrinos coming 
from the photomeson production are mainly expected. For 
the dashed curve where tj = \, pp neutrinos become impor- 
tant, as in the dissip ative photospheric scenario (iMurase 20081 : 
IWang & Dai 20091) and resulting in similar neutrino spectra, 
although our prediction comes from an independent motiva- 
tion so that a larger baryon loading is typically needed to have 
a visible ^ GeV component. For the dot-dashed curve, only 
pp neutrinos are relevant because inefficient acceleration is 
assumed, and we do not expect ^ PeV neutrinos because 
the proton maximum energy is not so high. In any case, for 
one burst at z ~ 1, the expected number of muon events in 
IceCube is ~^ 10" - 10"', which seems too small to detect. 
Hence, bright bursts with £™ > 10 erg or nearb y bursts at 
z '-^ 0.1 are typically necessary dPermer & Atovan 2003.) . al- 
though they would be rare. 

On the other hand, the cumulative neutrino background is 
of interest, since time and space coincidences are expected to 
be of use for analyses of GRB prompt emission. The cumu- 
lative neutrino background is calculated using the following 
formula. 
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where dN^/dE^ is the neutrino spectrum in the cosmological 
rest frame and R{z) is the GRB rate. For demonstration, we 
use the GRB 3 evolution but differences in evolution models 
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do not make significant changes (see lMurase 20071 and refer- 
ences therein). 

An example of the cumulative neutrino background ex- 
pected in the hadronic model is shown in Figure 8, using the 
parameters given in the caption of Figure 5. However, the ap- 
propriate normalization of the background flux is not so ob- 
vious, since £^ comes from the turbulent energy rather than 
the cosmic-ray energy in the re-acceleration/heating scenario, 
where S^^^/£^ can be small theoretically. Indeed, E-f^^/E^ 
should typically be smaller than the value inferred in Fig- 
ure 5 (e.g., Beniami ni et al. 201 ll). which can be achieved by 
changing sub-parameters (c.f. Figure 4). Therefore, on an ad 
hoc basis, the background flux is rescaled using £^ex instead 
of 8^. Since 8^^,^ is highly uncertain at present, for demon- 
stration purposes, we introduce the differential gamma-ray 
energy at 100 MeV, i^^ex, and rescale the background flux as- 
suming that it is a few % of the total gamma-ray energy. The 
muon neutrino flux is analytically estimated to be 
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where /^ is the correction due to the redshift evolution, 
and the characteristic neutrino energy is generally model- 
dependent (cf. Figure 7). One should also keep in mind 
the following limitations of the result. (1) We have used 
the one-zone neutrino spectrum, but contributions from larger 
emission radii would be dominated at high energies (e.g., 
iMurase & Nagataki 2006ah . (2) The background spectrum 
can be affected by i ndividual characteristics of the bursts (e.g., 
iBecker et al. 20061) . so that it would be valid for the purpose 
of estimating the flux level when typical parameters are cho- 
sen. In our case, if the extra hard component at > 10 MeV 
energies is less prominent for most of the bursts, the back- 
ground flux should also be less. Despite these caveats, the 
figure suggests that high-energy neutrino signals should be 
one of the important messengers of the hadronic model, and 
one may expect a few events per year in IceCube, provided the 
model assumptions are valid. Just for comparison, predictions 

in other scenarios are also shown in Figure 8. 

Recently, the IceCube collaboration jAbbasi et al. 201 lb 
has given upper limits on the neutrino background flux ^ 
(with a normalization procedure different from theoreti- 
cal work by Waxman & Bahcall 1997 and Murase & 
Nagataki 2006). Such limits are interesting in order to 
test th e GRB-UHECR proton h ypothesis for the prompt 
phase (I Waxman & Bahcall 1997h as the photomeson pro- 
duction efficiency of fp^ < 1 can be probed. Generally 
speaking, the efficiency is model-dependent. One typi- 
cally expects max[l,/p^] ^0.1-1 for emission radii of r ^ 
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lO'^^ cm (iMurase & Nagataki 2006ah . Especially in the neu- 
tron escape scenario where the neut rino flux is comparabl e 
to the escaping cosmic -ray flux (e.g.. lMannheimetal. 20001) . 
the limits are alread y stringent for steep proton spec- 
tra (lAhlers et al. 201 ll) . However, the detection of neutri- 
nos would typically be difficult if the ob served UHECRs are 
mainly heavy nuclei rath er than protons dMurase et al. 2008t 
IMurase & Beacom 20101) . Observations in the near future 
will give more stringent limits on fp^, based on the GRB- 
UHECR proton hypothesis. On the other hand, the hadronic 
model focused on in this work relies on the different moti- 
vation of explaining the observed GeV emission, where we 

' Note that time- and spa ce-coincidence betwe en neutrinos and gamma 
rays is assumed (but see, e.g.. lJacob & Prran 20071) 
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Fig. 8. — The cumulative neutrino backgrounds (for all flavors) from 
GRBs. Hadronic: calculated in this work for the demonstrative param- 
eter set used in Figure 5, with rescaled normalization of £hE7cx'^(0) = 
10 erg Mpc"-^ yr"'. Classical: originally predicted by Waxman & Bah- 
call (1997) and numerically calculated spectra (sets A and B) are taken 
from Murase & Nagataki (2006a), but with normalization of £'hecr^(0) = 
5 X 10^' erg Mpc"-^ yr"' . Photospheric: calculated in Murase (2008), but with 
normalization of gcR = ^cr/^7 = 1 and £™ = 10 erg. WB: the Waxman- 
Bahcall bound (Waxman & Bahcall 1 99^ shown as landmarks. MB: the ef- 
fective iron-survival bound tMurase & Beacom 201(1) shown for comparison. 



do not have to explain UHECRs by GRBs. Since it is un- 
certain how common and energetic the extra hard component 
at > 10 MeV is, it is not so obvious to derive quantitative 
constraints on the model. Nevertheless, although the back- 
ground flux seems lower than the current limits, our result 
implies that future neutrino observations would be crucial for 
the hadronic models provided that the extra hard component 
commonly carries 1 - 10 % of the MeV component. Limits at 
differential energies would especially be useful. 

5. SUMMARY AND DISCUSSION 

In this paper we studied the role of stochastic acceleration 
(or wave heating) on the GRB prompt emission, in the pres- 
ence of electromagnetic cascades. As an example, we employ 
the hadronic model to initiate cascade processes. We demon- 
strated that some of the current issues in the GRB prompt 
emission can be explained. (1) The low-energy photon index 
is expected to be a ^ 1 when electrons and positrons con- 
fined in the downstream region are accelerated stochastically 
or heated via turbulence, and the typical peak energy, £'''^ is 
stabilized by the balance between the stochastic acceleration 
(or heating) time and the cooling time. (2) It is possible to 
avoid the fast cooling problem because the charged leptons are 
slowly heated by turbulence. When the leptons are injected 
via the hadronic processes, we may not have a potential is- 
sue on injection of primary electrons in the shock acceleration 
mechanism. (3) An extra hard component, EF(E) oc /i"""^, 
which has often been observed above 1-100 MeV, can be ex- 
plained by cascade emission, such as in proton-induced cas- 
cades. The existence of the extra component is an important 
feature of the cascade initiated by high-energy injection well 
above MeV energies. 

Although we have considered one-zone cases, dissipa- 
tion will occur at various radii, which can affect the high- 
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energy spectrum. As shown in appendix B, in slow heating 
scenarios via stochastic acceleration, the high-energy pho- 
ton spectrum can be harder than in the one-zone case be- 
cause of the superposition of one-zone slow-heating spec- 
tra. As a result, the broadband spectral shape can be a bro- 
ken power law up to ~ GeV energies. Such multi-zone ef- 
fects are even more important above > 0.1 GeV. Smaller 
77 opacities at larger radii gene rally smear out the 77 pair- 
creation cutoff dAoiet al. 20101) . Note that, in our calcula- 
tion method (where spectra are obtained via iteration), the 
radiative transfer effect in the emitting slab, which leads to 
(1 — exp(— T^^))/T^-y rather than the exponential cutoff, is not 
included. Also, time-dependent and geometrical effects in- 
troduce additio nal compUcations, and the pair-creation break 
can be higher ( Granot et al. 20081) . Hence, we expect that the 
extra component extends to higher energies, although details 
are beyond the scope of this work. 

One of the observational properties found by Fermi is the 
delayed onset at GeV energies. It is not so obvious how to 
explain this in our model. However, there are several reasons 
that may lead to broader pulses at high energies. Especially, 
multi-zone effects seem important. As expected from the 
above, the cascade component as well as the slow heated com- 
ponent would be superposed by dissipation at various radii. 
Also, gamma rays from outer radii would be more transpar- 
ent, which could make a dominant contribution at GeV-TeV 
energies. In addition, in the hadronic models, neutron beams 
from the inner radius make further gamma rays via neutron 
decay and photomeson production with low-energy photons 
produced by dissipation at the outer radius. Leptons may also 
up-scatter the MeV emission from inner radii via the external 
IC process. In addition, in our model, emissions on timescales 
longer than fjur and fjyn are expected. Even after the dynamical 
time, the residual energy of electrons would be released by the 
external IC process and synchrotron radiation in some resid- 
ual magnetic field. In the hadronic models, further hadronic 
reactions may also happen as the escaping cosmic rays diffuse 
in the residual magnetic field (although they are adiabatically 
cooled). Hence, future refinements based on multi-zone cal- 
culations seem necessary. 

All hadronic models have typically several caveats or open 
issues. They often lead to the requirement of a high baryon 
loading, which may be problematic, although this is moti- 
vated mainly by the GRB-UHECR hypothesis. Another po- 
tential issue is the radiative efficiency. In our cases, we typi- 
cally obtain the reasonable radiative efficiency, L-y/Lth ^0.1, 
although the apparent radiative efficiency can be higher due 
to lower 7reacc, release of the cosmic rays, residual emission 
and subsequent internal collisions. Despite these concerns. 



there are interesting features. Importantly, the hadronic model 
seems testable by upcoming multi-messenger observations. 
The high-energy neutrino signal is one of the most impor- 
tant signatures, since the neutrino luminosity is expected to be 
comparable to the extra hard component. If this extra compo- 
nent is ubiquitous and it carries ^1-10% of the total gamma- 
ray energy, IceCube would detect this signal in multi-year ob- 
servations. It is also important to check whether the extra 
hard component is common among all GRBs. As of now, its 
detectability with Fermi is limited and the current analyses 
suggest that the GeV emission has less than ^ 10 % of the 
MeV emission. Therefore, in order to unravel the properties 
of the high-energy emission, deeper observations via detec- 
tions with Cherenkov detectors such as MAGIC and VERI- 
TAS should be important, even though detections from distant 
GRBs become difficult due to the attenuation by the extra- 
galactic background light. Although very-high-energy pho- 
tons from GRBs have not been fir mly detected so far, the fu- 
ture CTA (ICTA Consortium 20101) and HAWC ai-e anticipated 
to change this situation. 

We note that, although we have focused on the hadronic 
model as an example of high-energy injection, our results can 
be applied to other models, e.g., such as magnetic dissipation. 
In the slow heating scenario resulting from stochastic acceler- 
ation (or heating), it seems common to have hard low-energy 
photon indices, a^\ and /? > 2. This feature itself would not 
be changed by details of the high-energy injection, as long 
as slow heating can operate and the MeV peak comes from 
balance between acceleration/heating and cooling. Hence, as 
for the origin of stochastically accelerated (or heated) parti- 
cles, one may think of other high-energy injection processes 
and subsequent cascades. If magnetic dissipation such as re- 
connection accelerates electrons up to very high energies and 
subsequent cascades are developed, similar spectra would be 
obtained. 
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APPENDIX 
STOCHASTIC RE- ACCELERATION OF MESONS AND MUONS 

Here, we consider the possible effect of re-accelerated mesons on high-energy neutrino spectra. First, foxE^'^ < Ef,, the classical 
neutrino spectrum produced by shock accelerated proton s via the photomeson production is written as (.Waxman & Bahcall 1999t : 
iRachen & Meszaros 1998HWaxman 201lHMurase 2007h 



{Et/El'f~\EjElf 



{E,<El') 
{El' <E,< Et) 
(Et < E,) 



(Al) 



Here E^ is the low-e nergy break coming from E^' « Q.QSE^j or the characteristic energy where fpj(20E^) = 1 (' Asano 20051 : 
iMurase & Nagataki2006a) . Also, we have only assumed the pion synchrotron loss as the important cooling process, and E^ « 
0.25£'^ is the high-energy break coming from the pion synchrotron loss, though other processes such as adiabatic cooling could 
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also be relevant. From f^syn = fir = TttT^ (where Tjr is the lifetime of charged pions), the high-energy break energy of pions is 
estimated to be 



17-' ^-r / ^T^if^TrC^) 1QAD\/ -1/2 ,-1/2 T.2 /AOX 

E. « ry ^^^,^5g2,^ - 190 PeV .,,_, r,3L.,,53.5r2.. (A2) 

Note that one can make similar estimates for other particles such as kaon and muons, changing proper mass and lifetime. 

If there is re-acceleration by turbulence, f^.sta = fTr.syn gives the maximum energy of accelerated pions unless the energy fraction 
carried by mesons exceeds e^, similarly to the case of electrons. When one uses equation (9), the typical energy of pions is 
estimated to be 

6nmlc\eBf-'' 



eT'. 



armlB^ll^^ r/sm 



(A3) 



typ 



The situation depends on details of the plasma/MHD turbulence at relevant scales. As numerical examples, we have £7? 

120 PeV e^'/^i3^'Craoi33.5r2^7''7fm' for ? = 1-8, and eT c, 880 PeV e-,][',rni;^JJoL2^l,Ti,r^T for <? = 1-5. However, the 
re-acceleration of charged mesons is possible only when particles are injected to the re-acceleration process before their decay. 
The critical energy is determined by f^r = fjr.sta, and we find 




For 1 < ^ < 2, we expect two cases, E" < E^ < £'^'' and £'J?'' < E^ < E", and re-acceleration of mesons is relevant only in the 
former case at energies above E'^. For q < I, two cases, fij < E'^^'^ < E^ and E" < £'^'' < E-^, are possible, and mesons with 
energies below E" are re-accelerated in the former case. Usually such re-acceleration will not be important at large dissipation 
radii especially above the photospheric radius, since time scales such as fsta, fcooi and fdyn become longer for larger radii whereas 
t-jr is unaffected. If the meson re-acceleration happens, mesons are re-distributed so that the resulting meson spectrum will be 
altered. As long as particle escape is irrelevant, the neutrino spectrum can be hard, and one may expect 

eI(J),ocE^-HE,<eT). (A5) 

The typical neutrino energy can be very high energies of ^ 100 PeV, so that detectability of such very high-energy neutrinos may 
be enhanced by re-acceleration. 

SUPERPOSITION OF SLOW HEATING SPECTRA 

In the slow heating scenario with stochastic acceleration (or heating), whether the origin of relativistic electrons (and positrons) 
is leptonic or hadronic, one may expect a synchrotron component with a hard spectrum, 

EF(E)o<:E^-"(E<E^'), (Bl) 

where a ~ q/2 if the MHD turbulence is responsible for stochastic acceleration. Above £■*", the spectrum would typically have 
a cutoff feature or significant steepening {/3 ^ 2). 

So far, we have assumed the one-zone case, where dissipation occurs at the specific radius, r. But, realistically, dissipation 
continues at further radii, and we discuss its effect on the high-energy spectrum. For simplicity, let us here assume that Ly,LB oc 
Ldiss oc r"' (s > 0). For a given T ex r~^, assuming Ztur « / ex r/T, one obtains 

fsta CX B^-^-f^-nl' CX ^?-l-(9/2-l)(.'+2)«(2?-3)^2-,^ (g2) 

On the other hand, the (synchrotron/IC) cooling timescale has 
Then, for the synchrotron break energy, we have 

E'^'^OcTj^ ^^^(s9+29-4)/(3-9)+1-V2+2a:(1-29)/(3-?)^ ^-g^-j 

As a result, if the one-zone high-energy photon spectrum above £''' is steep enough, the multi-zone high-energy photon spectrum 
may have a harder one, 

EF(E) CX E 2(.„,+2,-4H.,-2)(3-<,)+441-2<,) . (B5) 

Especially if q ^ 2 and x ^ 0, one has EF(E) ex ^-'/Ci+i 5s) which may be responsible for the high-energy index. In this case, 
higher-energy gamma rays are emitted mainly from larger radii, so that those timescales can in principle be longer. Note that the 
above discussion can be applied to not only the shock dissipation models but also the magnetic dissipation models, and a broken 
power-law spectrum can be expected in either case. 
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